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Harmonically Related Bands 
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A method is proposed for the study of the solar corona, by observing a delayed correla- 
tion between rapid fluctuations of enhanced solar radio emission in harmonically related 
frequency bands. The correlation is expected in those types of emission that are produced 
by nonlinear plasma oscillations. The delay of the fundamental frequency with respect to 
the harmonic would be brought about by dispersive group retardation in the corona. The 
method appears to be most suitable for use with type II bursts, though it might also be 
applied to other types of nonthermal solar emission. 



A solar noise burst of typo II consists of a narrow 
band of intense radio emission, that drifts steadily 
downwards in frequency over a period of several 
minutes (Wild and Macready [1] and 1 Wild [2]). Often 
two such bands are observed simultaneously, with 
center frequencies roughly in the ratio 1:2. This 
fact suggests that the burst originates in nonlinear 
plasma oscillations in the solar corona (Wild, Murray, 
and Rowe [3]; Maxwell, Swamp, and Thompson [4]). 
The oscillations are thought to be excited by the 
motion through the corona of a bunch of particles 
that has been ejected violently from the photosphere. 
Usually the burst is preceded by a solar dare which, 
in these cases, may well be the original source of the 
fast particles. At each point along their course 
outward through the corona, the impact: of the 
particles excites oscillations at the local plasma 
frequency, which radiate the noise observed in the 
lower of the two bands. Because the oscillations are 
supposedly nonlinear, energy is radiated also at the 
second harmonic of the plasma frequency, and this 
corresponds to the upper hand. As the bunch ol 
particles continues outward, it. encounters progres- 
sively decreasing local plasma frequencies, and this 
circumstance explains the downward progression of 
the frequency of the burst. 

Now it is possible that the noise in these bursts 
may exhibit fairly rapid nonstatistical fluctuations. 
of duration 0.1 sec or less, due to irregularities either 
in the distribution of electrons in the corona or in the 
strength of the source of excitation for the oscilla- 
tions. If so, then the same fluctuations would be 
expected to occur on both the fundamental frequency 
and on the harmonic, since the two frequencies are 
generated simultaneously at the same point by the 
same process. However, the fluctuations would not 
be exactly simultaneous when observed at the earth, 
because, as the waves travel outward through the 
corona, the different frequencies would experience 
different amounts of group delay; the fundamental 
would be delayed more than the harmonic. If the 

i Figures in brackets indicate the literature references at the end of this paper. 



differential delay between the fluctuations on the 
two frequencies could be observed and measured, it 
might prove to be a source of useful information 
about the solar corona. 

To find out whether this differential delay | is 
Likely to be observable, its size has to be estimated. 
In the appendix, some calculations are made on the 
assumptions that the electrons in the corona} are 
distributed with spherical symmetry about! the 
center of the sun, and that their number density 
A r varies with distance from the center according 
to t he power law 



N=N Ii- b , 



(1) 



where /Y is the nominal density at the surface of 
the photosphere, and the distance R is measured in 
so In i- radii. The source of excitation for the burst is 
assumed to be traveling radially outwards from 
the sun, along the line directly towards the earth. 
In calculating the group delays, no account was 
taken of the possibility that magnetic fields might 
be present. On these assumptions, the differential 
delay between a fluctuation on the fundamental 
frequency/ and the same fluctuation on the second 
harmonic (2/) is 



AT(f)=Kf- t,b , 



(2) 



where the factor of proportionality K depends on 
the values of N and b. These parameters were 
estimated by fitting the distribution (1) to the 
densities given by Van de Hulst [5], which were 
obtained experimentally by measurement of the 
brightness and polarization of coronal light at dif- 
ferent angular distances from the sun. The densities 
used were those that apply at sunspot maximum, 
and the fit was made at the level where the density 
corresponds to a plasma frequency of about 100 
Me; this process gave the values iVo=4X10 8 cm -3 , 
and 6 = 9. Thence i£~l, when AT is measured in 
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seconds, and / in megacycles. Equation 2 then 
yields the following estimates of the delay at various 
frequencies of the fundamental: 



/ 



AT 



80 Mc 


0.40 sec 


100 Mc 


.38 sec 


120 Mc 


.36 sec 



These are lower estimates because, above the 
100 Mc level, the densities given by (1) are rather 
less than those of Van de Hulst's distribution; yet, if 
that distribution is correct, then they are not likely 
to be excessively low, as most of the differential 
delay of the waves is produced just above their level 
of origin. 2 

It appears that a differential delay of about 0.4 
sec should be expected, and that, for the assumed 
distribution, it should vary only slowly with fre- 
quency. A delay of this order should be observable 
readily, provided that the noise itself contains 
fluctuations of much shorter duration. If observed 
on just one harmonically related pair of frequencies, 
the delay would give a weighted measure of the total 
electron content of the corona above the level of 
origin of the waves. If the delay could be measured 
at several such pairs of frequencies, and the same 
assumptions were made as in the calculations above, 
then the observed variation with frequency would 
give the value of the index b in the power law for 
the distribution of electrons. More generally, it 
should be possible to infer the distribution of electrons 
in height with no prior assumptions about its form, 
by a variant of the "lamination" methods that have 
been developed for analyzing the group delays 
observed in ionospheric sounding (Thomas [8]). 
This information, when combined with the observed 
frequency/time relation for the burst, would show 
how the height of the source of excitation in the 
corona increased with time. At present, the radial 
movement of the source can be traced only by 
assuming some model for the distribution. 

In this kind of analysis, height would be measured 
with respect to the level of origin of the highest 
fundamental frequency that was observed. The 
height of this particular level with respect to the 
photosphere would have to be estimated in some 
other way. For instance, if the radio burst was 
preceded by a flare, then the time interval might be 
measured between the occurrence of the flare and 
the onset of the burst; this information, when com- 
bined with the calculated velocity of movement of 
the source of excitation, would give the height. 

Such an analysis would also be complicated by 
the fact that the source would not, in general, be 
moving directly towards the earth; nor would it 



2 However, some recent studies indicate that Van de Hulst's figures may be 
too low. Newkirk [6], on the basis of further optical measurements, concludes 
that even on the quiet disk the figures should be multiplied by a factor of 2.2, 
while over active regions the factor should be 4.5. Maxwell [7], on evidence from 
the radio bursts themselves, has suggested that the multiplying factor for active 
regions may be as high as nine. If such factors indeed apply, then the estimates 
of the differential delay A7 1 should be multiplied by these same factors raised to 
the power of 1/6 (see eq 11 in the appendix); for the given value of b, the increase 
in the estimates is less than 25 percent. 



necessarily, even, be moving radially outwards from 
the sun (Wild, Murray, and Rowe' [3]). However, 
it is worth noting that any motion of the source 
transverse to the line of sight can be observed 
directly by other means (Payne-Scott and Little 
[7]). If such observations were combined with those 
of the differential delay, it might then be possible to 
infer the distribution on the assumption of spherical 
symmetry. This assumption itself may seem naive, 
considering that the corona is known to have an 
irregular structure, but nonetheless it may still be 
useful as a basis for deriving a first approximation 
to the distribution. 

The first step in implementing these proposals 
would be to find out whether the noise in the type II 
bursts does exhibit n on statistical fluctuations that 
are sufficiently rapid to be of use for measuring the 
differential delay. This question cannot be an- 
swered from the results of previous observations, 
since all these were made using receivers with too 
large band widths. In detecting a signal that varies 
in frequency, the best bandwidth to use is roughly 
equal to the square root of the rate of change of fre- 
quency (Storey and Grierson [10]). Judging from 
the spectral characteristics of the type II bursts, 
the best frequency range in which to observe them 
on the fundamental would be 80 to 120 Mc, and in 
this range their rate of change of frequency is 
usually about 1 Mc 2 . Hence the RF bandwidth of 
the receiver should be about 1 kc. With so narrow 
a bandwidth, a fairly large antenna would be required 
to provide adequate sensitivity. 

If rapid fluctuations were observed, the next step 
would be to search for the anticipated correlation 
between the fluctuations in harmonically related 
bands. The antenna would be connected to a pair 
of receivers, one tuned to a fixed frequency in band 
80 to 120 Mc, and the other to twice this frequency; 
the high-frequency receiver should have an RF band- 
width double that of the low-frequency receiver (the 
criterion of rate of change of frequency would require 
the bandwidth at the harmonic to be ^2 times that 
at the fundamental, but this requirement is not at 
all critical, so it seems better to arrange that the 
limits of the two bands are related harmonically). 
The outputs from the detectors of the two receivers 
would be recorded on magnetic tape, together with 
timing signals. If a particular section of tape was 
identified as the record of a type II burst, then this 
record would be played back at reduced speed, the 
information on the two data tracks converted to 
digital form, and the cross-correlogram of the two 
sets of fluctuations evaluated in a digital computer. 
In this way a search would be made for a correlation 
between the fluctuations observed in the two bands, 
and for a possible systematic time displacement that 
would give the best correlation . If this search were 
to prove successful, then the measurements would 
be extended to cover more than one pair of fre- 
quencies. There would be advantage to making 
these measurements in conjunction with recordings 
on a panoramic radio spectrograph, which could be 
used to identify the type II bursts. 
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Possibly the same techniques could be used to 
study other types of enhanced solar radiation. Thus 
the noise bursts of type III also exhibit harmonic 
struct n iv (Wild, Murray, and Rowe [3]), but the rate 
of change of frequency is much faster, of the order of 
50 Mc 2 ; hence the receiver would need to have a 
rather larger RV bandwidth, about 7 kc, for recording 
these bursts. A further possibility is the study of 
noise storms. Though these disturbances are not 
localized in frequency at any one time, as are the 
bursts, but to the contrary are spread over a wide 
frequency band (Wild and Macready [1]), neverthe- 
less they also may be generated by the same process 
involving nonlinear plasma oscillations. The differ- 
ence between the two kinds of disturbance may lie 
only in the nature of the excitation; a continuous 
stream of particles for the storm, as against a discrete 
bunch for the burst. If so, then there would be the 
same grounds for expecting a delayed correlation 
between fluctuations of noise power in harmonically 
related bands; in this case, however, there is no way 
of telling in advance what RF bandwidth would be 
best for observing the correlation. I am indebted to 
R. S. Lawrence for this suggestion. 

In this work, there is one possible source of diffi- 
culty that needs to be recognized. It is the presence 
of irregularities in the corona, which may cause the 
radiation to be propagated from the source to the 
earth by more than one path. Multipath propaga- 
tion could give rise to rapid and severe fading that 
would not be correlated in harmonically related 
frequency bands, and any correlated fluctuations 
might thus be masked. A similar difficulty has been 
met with already in an attempt to measure the 
electron density at great heights above the earth by 
observing the dispersion of radar echoes from the 
moon (Eshleman, Gallagher, and Barthle [11]). 
However, it is hard to judge whether this factor will 
present a serious obstacle to the proposed studies. 



In writing this paper, I have been helped con- 
siderably by advice from Dr. Constance S. Warwick 
of this Laboratory, and from Dr. G. Newkirk of the 
High Altitude Observatory of the University of 
Colorado. 
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Delay 



In calculating the differential delay, the following 
symbols will be used to refer to points, first, at the 
surface of the photosphere, and second, at some level 
in the corona: 

Photosphere Corona 

Distance from center of sun r r 

Local electron density N N 

Local plasma frequency F j v 

The electron density in the corona will be assumed 
to depend on r only, according to the power law 



N=N R~ b 



(1) 



where 



R=r/r 



(3) 



The plasma frequency /„ is proportional to the 
square root of N, so its distribution is 



fp=foR-» 2 . 



(4) 



By hypothesis, the oscillations that produce the 
noise emissions take place at the local plasma 
frequency; emissions of frequency / are produced 
at the level, r(f) say, where f p =f. The difference in 
the group delay experienced by the fundamental of 
frequency / and the harmonic of frequency 2/, in 
traveling from this level out through the corona to 
the earth, is 



AT(/) = l/cf 00 {//(f)-,/ (2/) dr 

Jr(f) 



(5) 



where c is the speed of light, and // is the group 
refractive index. If magnetic fields are absent, as 
will be assumed, then // depends only on the relative 
values of the wave frequency and local plasma 
frequency, according to the law 



where 



p'(f) = (l-X)-l 

x=ujf)* 



(6) 
(7) 



Along the path of the waves, the dimensionless 
quantity X varies from a value of unity at the source 
to zero at infinity. It is convenient to adopt X as 
the variable of integration: 

AT(/) = l/cf°{M / (/)-M / (2/)}(^)dZ 

= l/eJ o 1 {M'(/)-M , (2/)}(-^)rfX. (8) 

The components of the integrand art 4 obtained as 
follows: first, from (6), 

^(f)-/(2f) = (l-X)-i-(l-X/4)-l (9) 

Also, from (7), (4), and (3), 

dr 



-J=^(/o/i) 2/6 i- (1+1/J) - 

Hence the differential delay is 

AT(f)=Kt 2/b 



whei 



be 



(10) 

(2) 
(11) 



and 



1= i l {(l-X)-i-(l-Xl4)-i}X- a+l "'\lX. (12) 
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Now the factor K must be evaluated. The con- 
stants are c=3X 10 5 kmps, and r =7X 10 5 km, while, 
from the data of Van de Hulst [5], 6^9 and / — 180 
Mc.D The integral / can be evaluated numerically, 
though with some difficulty, since one or other of the 
factors in the integrand becomes infinite at each limit. 
The result is that, when / is measured in megacvcles ; 
#=1.(1). 
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